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In an earlier work [P. J. Bardroff and S. Stenholm, submitted to Phys. Rev. Lett.] , we have derived 
a fully quantum mechanical description of excess noise in strongly damped lasers. This theory is 
used here to derive the corresponding quantum Langevin equations. Taking the semi-classical limit 
of these we are able to regain the starting point of Siegman's treatment of excess noise [Phys. Rev. 
A 39, 1253 (1989)]. Our results essentially constitute a quantum derivation of his theory and allow 
some generalizations. 

■q^- . PACS numbers: 42.50.Lc, 42.55.Ah 
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I According to the fluctuation-dissipation theorem, the attenuation or amplification of a signal always adds noise. 
In optical amplifiers, this fact is usually phrased as "one noise photon" added to the signal from the spontaneous 
emission processes in the reservoir. This assumption about the noise gives rise to the phase diffusion responsible for 
the Schawlow-Townes linewidth of lasers. However, this is not true generally. In particular cases, the noise can 
exceed the intensity of one photon by the so-called excess- noise factor or Petermann ii'-factor [Q. 

Experimentally this phenomenon was first confirmed in a laser cavity with large output coupling leading to an 
enhancement of a few times Later, even a factor of a few hundreds was achieved for solid state lasers and 
• gas lasers . Also experiments with a coupling of the polarizations and an inserted small aperture ||^ have 
] demonstrated large excess noise. A recent experiment has shown that excess noise can be colored due to saturation 
O • effects 0. 

2^ ' After the prediction of excess noise by Petermann for the case of gain-guided semiconductor lasers Q, the first 
more general theory of excess noise was given by Siegman using a semi-classical description [0. Until recently, only 
a few simple systems have been discussed from a quantum mechanical point of view |12[| . 

In a previous paper fl^ we introduced a master equation describing a multi-mode field interacting with a reservoir 
' describing the general linear amplifier or attenuator in a strictly quantum mechanical formulation. We find that under 
^ ' certain conditions, the reservoirs create couplings between the undamped modes of the system. Such dissipative 
^ couplings lead to a non-Hermitian eigenvalue problem, which introduces non-orthogonal quasi modes in a natural 
manner. The amplitudes of these modes are then found to display the expected excess noise, which we here ascribe 
^ ^ to the reservoir-induced mode-mode coupling. 
•'H . In this paper we derive the quantum Langevin formalism following from our theory in ReL Wher eas the 

r^J ' dynamic variable of the master equation is the quantum state, the Langevin equations are for the field operators. 
■ This allows us a direct comparison of our approach with the well-known semi-classical treatment introduced by 
Siegman . As this has provided the physical understanding and the mathematical expressions for the excess noise, 
we are pleased that we can essentially derive his starting equations from our fully quantum mechanical treatment. 
We are also able to generalize the semi-classical analysis of excess noise to cases beyond the paraxial approximation. 



II. MASTER EQUATION 

In this section we briefiy review the results of the quantum derivation of the excess-noise factor based upon the 
master equation. We use orthonormal real mode functions Un{x) of the electromagnetic field with frequency a;„ which 
fulfill the boundary conditions for the given configuration in the whole "universe" and satisfy the orthonormality 
relation 

J SxUn{x)Um{x) ^ 5nm, (1) 

where V is the volume of the whole space. Note that the mode function Un{x) is a vector including the polarization 
orientation and that we choose them to be real for convenience. The electric field operator then reads 
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E{x) = '^enUn{x) (ctn + al) , (2) 

n 

where a„ and ajj are the usual creation and annihilation operators of the field excitations and the so-called vacuum 
field amplitude is 
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We start from the multi-mode master equation 
d 



(3) 



+ '^^^ni,n{'2dnp{t)ai^- al^anp{t) ~ p{t)a'l^an} - i^uj„[alan,p{t)]. (4) 

n,m n 

with the two symmetric matrices Tm.n and Lm.n given by 
and 
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The former describes losses and the latter amplification due to the interaction with the reservoirs. 

The two reservoirs for amplification and attenuation are assumed to consist of two-level atoms injected in the upper 
or lower state, respectively |2]. They are completely characterized by the position dependent injection rates ri(x) and 
rr (a;) , the interaction time r of the individual atoms with the field and the orientation of the atomic dipole moment 
d. In principle, the dipole orientation could be different for damping and attenuation and it may depend on position. 
This treatment of the damping can describe spatially localized absorption due to an inserted aperture or due to a 
detector placed outside the cavity. Assuming a perfect absorber (or detector) surrounding our cavity, the reservoir 
can also model the damping due to output coupling. Taking the limit of a infinitely large "universe" (V — > oo), and 
hence using a continuum of modes, would be another way of including losses due to output coupling in our model as 
shown in Ref . @ . 

Because of the interaction through the reservoir, the time evolution of the mean values 

m 

exhibits coupling between different modes. The definition of the quasi modes operator A follows from imposing the 
condition 



d 
di 



{A}^{^{X-^)-^n]{A}, (7) 



where n is the frequency, A is the amplification rate and 7 is the attenuation rate. For later convenience, we split the 
net-amplification rate (A — 7) into the two separate contributions A and 7. Note that O, A and 7 are real. We write 
this mode operator in terms of the free field mode operators as 

f i = ^ £„c„a„ (8) 

n 

with the expansion coefficients c„. This transformation includes the vacuum-field amplitudes e„ and we define 
£ — ^J~^^, because then the classical field amplitudes £n{dn) obey the same transformation as the operators. 
Inserting Eq. (^ into (^) we get an eigenvalue equation 



i) - i5n,m^n \ J^C„ = | ^ (A - 7) - if] [■ (9) 
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for the non-Hermitian matrix {■^{Lm,n — ^m.n) — *^n,m'^n}f^- Here Cn*^^ is the right eigenvector; the corresponding 

left eigenvector is s'^ch '' The superscript v distinguishes the different eigenvectors. The detailed properties of 

the quasi modes are summarized in the Appendix. 

We can now calculate the noise of the quadrature operator 



(10) 



with the definition of the quasi- mode function U^{x) given by Eq. ( [A7| ) in the Appendix. Taking the noise averaged 
over position and comparing to the usual single mode master equation with the same frequency 17^, damping rate 7^ 
and amplification rate Aj/ we find an enhancement by the factor 
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Ji^) |2 



p2 ^(•^y 



(11) 



m,n and Vm,n, defined 



for the noise added by the reservoir; cf. Ref. ||13|. The excess noise is large when the matrices L 

in Eqs. (^), have large off-diagonal terms and when they are not identical. The former follows when the injection 
rates rL{x) and rr(a;) are not spatially constant whereas the latter when damping and amplification are spatially 
separated. 



III. QUANTUM LANGEVIN EQUATION 

Following the usual treatment we replace the time evolution described by the master equation (jj) by an 
equivalent quantum Langevin equation. This contains non-commuting noise forces which are designed such as to give 
the same moments as those derived from the master equation. The quantum Langevin equation is written 

^a(*)n = ^ X! ^ ^m-n) a,„i{t) ^ «t^na„(t) + fn{t), (12) 

where the Langevin noise sources fn{t) obey the correlation relations 

{Ut)fl{t'))^2{D^^^,J6{t^t'), (13a) 

{f2{t)Ut'))^2{D^,^^J6{t-t') (13b) 

and 

iUt)) = {Ut)ut')) = {fi{t)fi{t')) = 0. (14) 

It then follows from Eq. ( p^ ) that the expectation values obey Eq. (H). The diffusion coefficients (-D. -t ) and (Z).t ^ ) 

have to be determined to give the correct noise correlations (aJjOm) and {ama\)- We compare the time evolution of 
the noise correlations derived from the master equation, given by Eqs. (A12), to the one derived from the Langevin 
equation (n2h to find the relations 



{fi{t)am{t) + ai(t)fra{t)) = Tm,n (15a) 

and 

{fn{t)al{t) + an{t)fL{t)) = L™,„. (15b) 

With the help of the Einstein relations 

2(^a„at ) = |(a™(04(0) - («™(0(|4(i) " flm - ((|«™(i) - UmlM), (16a) 

^(^alaj = |(4(i)«m(0) ^ ((^^(O " /I(i))«™(i)) (4(t)(|«™(i) " /™(t))), (16b) 
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and Eqs. (|T|), we determine the diffusion coefRcients to be 

2(^a„at > = W + 4 W/m W) = r„,„ (17a) 

and 

2(^a1;a„) = ifrM&Ut) + W) = Lm,n. (17b) 

These relations clearly show the mode correlations due to the reservoir. 

In the following, we derive a wave equation for the propagation of the electric field operator including amplification, 
damping and the corresponding noise source. Starting from Eqs. (|l|) and (|) we can find the exact equation 

f 9^ 1 ( d ^ 

< - c^V^ \ E{x,t) - ^ e„w„(a;)| Y^{Lk.n - Tk,n)-^ak + h.c.j = 

^ ■' n k 

e„-u„(a;)| - ^ ^(Lfc,„ - Tk,n){Li,k - ^i,k)ai 

n k.l 

k 

k 

with the mode functions Uri{x) fulfiUing the Helmholtz equation 

{c^V"" + ujl)ur,{x) = Q (19) 

together with the appropriate boundary conditions. At this point we introduce a number of approximations based 
on the assumption that the average oscillation frequency Q of the electric field is much higher than the decay or 
amplification rates, e.g. |Ln,m| or |r„^m|. This is well justified in the optical regime where the former is at least six 
orders of magnitude larger than the latter. The spectral width Aw of the relevant frequencies uon is assumed to be 
of the order of the decay or amplification rate. To be more specific, we will neglect terms of the order 0(A,y/(D)^, 
©(Ai/Aw/w^) and ^^(Aw/cD)^ or smaller, and we assume 0{Xv) ~ C(7^)- 

On the LHS of E g. jl8| ) the two terms ■^E{x, t) and c^V^£^(x, t) are of the order 0((D)^. Since the remaining term 
on the LHS of Eq. (18) is proportional to the damping and amplification rate, we can approximate the frequencies by 
the mean frequency iD. Hence inserting the definitions of L^^rn and r„ „i, Eqs. (|^), and of e„, Eq. (H), the remaining 
term on the LHS of Eq. (^ yields 

^e„u„(x)^(Lfc,„ -Tk,n)^a,k + h.c. = 

n k 

^EnUnix) ^ ^£n£fc— / d^x' (r/, (x') - rrix'))[un{x')d][uk{x')d]—dk + h.c. = 

n k 

^Y^^nM^)^ j d^x'{rL{x')^rr{x'))[ur,{x')d]d^j^E{x',t)^ 

n 

2 - p o f) 

^ ' d^x'{rL{x')~rT{x'))5T{x-x')d®d^—E{x',t) = {Rl{x) ~ RT{x)) — E{x,t). (20) 



2eo7iy ' ' ' / - dt ' ' ' ' ' "dt 

Here the matrices L„^m and r„.m occur in their position representations 

^ ^ Unix) ® ul^{x')Lr,^m ~ ^7^*^^ (^)'^T ~ ^')'^ ® d^ = Rl{x)6t{x - x') (21a) 

n,m 

and 

X! '^"(^) ® Wm(a;')rn,m ~ ^^rrix)STix - x')d ® d^ = Rt{x)6t{x - x). (21b) 



n.m 
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Note that Rl [x) , Rr (x) , d® (f^ and the transverse (5- function 

5T{x-x') = ^^Un{x)®ul{x') (22) 

n 

are tensors. We can neglect the terms on the RHS of Eq. ( p^ ) containing the field operator since they are of the order 
0(Ai/)^ or ©(Ai/Aw), respectively. For the noise we only take terms of lowest order. Therefore we may neglect the 
first of the noise terms in Eq. (jl^) and we approximate « —iCjfn- Introducing the position representation 

=^£„u„(a:)/„(t) (23) 

n 

of the noise source we find 

^ - c^V^ - (RLix) - Rrix))-^'^ E{x,t) = -2iujf{x,t) + h.c. (24) 
The correlations of the noise operators are 

(/(x, t)f\x', t')) - ^Rr{x)6T{x - x')5{t - t') (25a) 



and 



Consequently the total noise on the RHS of Eq. (E4[) obeys 



{f\x, t)f{x\ t')) = ^Rl{x)5t{x - x')5{t - t'). (25b) 



{{-2iLof{x, t) + h.c.f) = {Rl{x) + Rt{x))5t{x - x')5{t - t'). (26) 

As expected, the effects of amplification and damping add for the noise whereas they subtract for the amplification. 



IV. SEMI-CLASSICAL TREATMENT 

Starting from Eq. (|24) we can now perform the transition to the semi-classical treatment replacing operators with 



c-numbers. The solution of Eq. (|2J) can conveniently be written using the positive frequency part E'^^^ of the 
electromagnetic field. The real part is the electric field and the imaginary part relates to the magnetic field. With 
the help of the Green function 

GW(x,x',i) =^{/,(a;)J7,(a;')e3(^--^")*-*f^-* (27) 

and the accumulated noise 

t 

N^+\x,t) - -2iuj J dt';^e^(A.-7.)(t-*')-^a„(t-t')[/^(2.)^ J d'x'U,{x')f{x\t'), (28) 







we find the field to be given by 

E^+\x,t)^^ j d^x'G'^+\x,x',t)E^+\x\Q) + N^+\x,t) (29) 
starting from the initial field £'^+^(x',t = 0). Within the approximations made, the quasi-mode functions Ui,{x) and 



U^{x) and their eigenvalues Ajy, 71, and f2i, are the same as defined using the master equation, Eqs. (A7) and ( |AS| ). 
When we now calculate the variance of the electric field E{x, t) averaged over position and compare with damping and 
amplification processes described by the usual single mode master equation, we recover the same _fC- factor, Eq. (pj]), 
as before. We find for the noise term 
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liK + - 7i. - 7^) - - 



V 



V 



X- / £xU^{x)U;{x)- / d^xU^{x){RL{x) + Rrix))U;{x). 



(30) 



Considering only one quasi mode, the noise simplifies to 



1 / d'x{{Nl+\x,t)+Nl-\x,t))^) 



huj^ exp{(A^ - 7^.)t} - 1 



(A^ + "f^)K^ 



£0 K - Iv 

with the enhancement factor in the commonly used form 

/ d^xU^{x)U;{x) J d^xUAx)U*{x) 



1/ d^xU^{x)U^{x) 



(31) 



(32) 



Note that our choice of the normalization for the quasi-mode functions is given as in Eq. ( [Aq ) . We have shown in 
Ref. [ID that Eq. (|3|) agrees with Eq. (|ll|) up to the order 0(Aw/u;)2. 

Siegman used an equation analogous to Eq. (|^ as the starting point for his derivation of the excess-noise 
factor. However, there are two interesting differences in the details of the noise source correlations, Eq. (p6|). 

The first difference in Ref. is that the spatial transverse (5-function is replaced by a usual (5-function and that 
the temporal 5-function is replaced by the Hertzian bandwidth Aw/(27r) of the reservoir. The latter circumstance 
is explained by using the Fourier representation of our noise correlation in Eq. (p6|); this leads to the same equation 
with 5{t — t') replaced by 5{u) — uo')/2n. We then integrate with respect to lu and lo' over the frequency bandwidth 
Aw to obtain 



dwdijj' 
27r 



5(w - w')e-*('^*-'^'*') + c.c. 



^ sin(Aa^(t-OMe--(t-*') 

2Tl{t - t') 



"2^ 



(33) 



for |f — t'l < (Aw) ^. This approximation is reasonable when the mean frequency w of the noise is much larger then 
the bandwidth Aw. 

The second difference in Ref. Q is that {Rl{x) + Rrix)) is replaced by 2 (75;^^) {Rl - Rr) = 2i?L with 
spatially constant Rl and Rr- This simplification is justified only when averaging over the whole volume V and when 
amplification and damping are balanced. For the derivation of the X-factor which involves an average over position, 
this is a valid replacement. However, one has to be aware of the subtlety that only non-constant Rl{x) and Rr{x) 
with Rl{x) 7^ Rr{x) lead to non-orthogonal quasi modes and hence can give K > 1. 

It is interesting to note that within the paraxial approximation, we obtain an equation analogous to the position 
representation of Eq. (jl^) — the starting point of our semi-classical analysis. Making the ansatz 

(34) 



E^+'>{x,t) = e'"(^/=-*)i?(+)(x) 
with E^^'>{x) slowly varying with respect to the longitudinal coordinate z, we get from Eq. (p^) 



c|-^(+)(x) = lliRLix) - Rrix)) + \ E^+Hx) + J{x) 



dz 



(35) 



where f{x,t) « e 



f{x). The time derivative d/dt of Eq. ( p^ ) is replaced by the derivative with respect the 
longitudinal coordinate cd/dz which is equivalent in a frame moving with the electromagnetic wave. The frequency 
part of Eq. ([l^) is replaced by the transverse Laplacian. Frequently, Eq. ( p5| ) is solved with mode functions of the 
transverse Laplace equation, depending only parametrically on the longitudinal coordinate. This distinction between 
longitudinal and transverse coordinates leads to a factorization of the ii'-factor into a longitudinal and a transverse 
part. 
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V. DISCUSSION 



In an earlier paper |13|| , we derived the master equation for a set of modes coupled to amplifying and attenuating 
reservoirs. This introduces couplings between the undamped modes of the total "universe" and leads directly to the 
introduction of quasi modes, which are found to exhibit the excess noise described originally by Petermann Q. 

Our treatment has been carried out only in the linear regime so far. This describes an amplifier or an attenuator, 
where the treatment is most straightforward and the results display the most transparent physical insight. How- 
ever, the saturation in an operating laser will need to be considered, and we are for the moment carrying out such 
calculations, which show the influence of the excess noise in the strong field situation. 

The best physical picture of this noise was provided by Siegman , who also supplied the quasi-mode expression 
for the excess- noise factor. This has then been used successfully to describe the experimental findings ||^-|^. In [T^ ] 
we showed that our quantum mechanical approach naturally provides an expression which is essentially identical with 
Siegman's results. 

Siegman, however, utilized a semi-classical Langevin approach, where the noise forces were added ad hoc to the 
classical equations for the amplitudes; the noise forces were then supplied with properly chosen correlation properties, 
which was shown to imply the presence of excess noise. Because this approach has been found to give both a physically 
attractive and theoretically justified description of the situation, we find it interesting to connect that treatment to 
our quantum approach in some detail. 

In this paper we derive the quantum Langevin equations following from our general master equation. Here we utilize 
techniques known from quantum noise theory, and obtain results that can be directly compared with the treatment 
of Siegman's, when the semi-classical limit is taken. Except for some minor differences, our resulting equations are 
identical with those used by Siegman. We thus claim that we have justified his formulation of the problem from a more 
fundamental quantum mechanical point of view. The differences found are either based on natural approximations or 
obvious qualifications of the results as e.g. the introduction of the transverse delta function in the noise correlations. 
In addition, we have been able to generalize the theory to situations outside the paraxial approximation. 

The results of our treatment, however, have bearings beyond the problem of excess noise in highly lossy cavities. The 
approach is quite general, and in addition to the Markov approximation we only need the rotating wave approximation 
for the interaction with the reservoirs. The master equation is then derived from first principles, and the nonorthogonal 
quasi modes emerge in a natural manner. The theory is fully general and may well be applicable to other high loss 
physical systems as well. For the moment we know of no observation that would show the equivalent of the laser 
excess noise, but novel situations may soon turn up. The lively research activity in quantum information processing, 
atom optics and novel measurement situations may provide potential applications of the present theory. 

The physics of our approach resides in the coupling of the undamped modes through the reservoirs. In such a 
situation, the only essential assumption in our derivation is the Markovian approximation. In highly damped systems, 
this may not necessarily hold, and the introduction of memory effects in our theory has not been considered so far. 
Some features are, however, expected to survive, but also unexpected complications may appear. These questions 
remain to be investigated. 



In this Appendix we recall from Ref. [|l3| those properties of the quasi modes which are relevant for the derivation 
of the excess noise. The only properties of the left and right eigenvectors of non-Hermitian matrices which we need for 
our analysis are their mutual orthogonality and completeness [ |l6| : The eigenvectors fulfill the orthogonality condition 




APPENDIX A: PROPERTIES OF THE QUASI MODES 





(Al) 



n 



n 



and the completeness relation 




(A2) 



with Y.n' ^I'^'n' ^ 0- Therefore, we can uniquely define the set of quasi-mode operators as 
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with the vacuum field amplitude 



A - — Vc^^^e a 
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The inverse transformation is 



Consequently the positive frequency part of the electric field operator is given by 



The quasi-mode eigenfunctions 



2 (^) 



UAx) = V " " , . Unix) 



satisiy an orthogonality relation 

V 

with their adjoint quasi-mode functions 

U.{x) = Y,c(:^Un{a 



The properties 



A 



and 



E=-2|„('^)|2 h \^ p2|„('^)|2' 

If 



can be obtained from the real and imaginary parts of Eq. after taking the scalar product with the vector 
From the master equation (||) follows the time evolution of the noise correlations 



k k 



and 

+iiuJn - UJrn){amai) + T 



For the quasi-mode operators, we find the correlations to be 



d_ 
di 



{AlA^) = |i(A, + - 7. - 7^) + - (44) + ^ ^im,n£n£r„c(r)*c(^) (Al3a) 



and 



^(44) = |i(A, + - 7. - 7/.) + ^{^. - (44) + ^ E ^m,nSnemci:^*c<;^l (A13b) 
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